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Abstract. Several properties of massive stars with large 
effects of rotation and radiation are studied. For stars with 
shellular rotation, i.e. stars with a co nstant angular veloc- 
ity f2 on horizontal surfaces (cf. Zahn 1992 ), we show that 
the equation of stellar surface has no significant departures 
with respect to the Roche model; high radiation pressure 
does not modify this property. Also, we note that con- 
trarily to some current expressions, the correct Eddington 
factors r in a rotating star explicitely depend on rotation. 
As a consequence, the maximum possible stellar luminos- 
ity is reduced by rotation. 

We show that there are 2 roots for the equation giving 
the rotational velocities at break-up: 1) The usual solu- 
tion, which is shown to apply when the Eddington ratio 
r of the star is smaller than formally 0.639. 2) Above this 
value of r, there is a second root, inferior to the first one, 
for the break-up velocity. This second solution tends to 
zero, when T tends towards 1. This second root results 
from the interplay of radiation and rotation, and in par- 
ticular from the reduction by rotation of the effective mass 
in the local Eddington factor. The analysis made here 
should hopefully clarify a recent debate between Langer 
( |l99"7| , |l99l| ) and Glatzel (|l99§). 

The expression for the global mass loss-rates is a func- 
tion of both Q and T, and this may give raise to ex- 
treme mass loss-rates (fir-limit). In particular, for O- 
type stars, LBV stars, supergiants and Wolf-Rayet stars, 
even slow rotation may dramatically enhance the mass loss 
rates. Numerical examples in the range of 9 to 120 M Q at 
various T e ff are given. 

Mass loss from rotating stars is anisotropic. Polar ejec- 
tion is favoured by the higher T e g at the polar caps (g e fi- 
effect), while the ejection of an equatorial ring is favoured 
by the opacity effect (^-effect), if the opacity grows fastly 
for decreasing T e g. 
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1. Introduction 

Recent models of stellar evolution with rotation (Meynet 
and Maeder 2000) have shown that rotation heavily mod- 
ifies all the model outputs for massive stars. In the course 
of the above mentioned work, it was realized that several 
basic points in the stellar physics need to be further clar- 
ified, since they may have some important consequences 
on the evolution. These points concern in particular the 
correct expression of the break-up velocities and the de- 
pendence of the mass loss rates M on the observed rota- 
tion velocities v. These problems are of great concern for 
the most luminous stars close to the Eddington limit, like 
OB stars, supergiants, LBV and WR stars. 

There is an interesting debate in recent litterature 
about what is the correct expression for the critical ve- 
locity and what is the dependence of the M-rates on the 
rotation velocities v. The critical rotation velocity of a star 
is often written as i^ rit = 2^(1 — T), where T = L/L Edd is 
the ratio of the stellar luminosity to the Eddington lumi- 
nosity. With this expression, Langer ( [1997 , 1998| ) suggests 
that "no matter the rotation rate may be, it (the star) 
will arrive at critical rotation well before T = 1 is actu- 
ally reached". Consequently, Langer introduces the con- 
cept that the stars generally reach the break-up limit, i.e. 
the f2-limit, earlier in evolution than the T-limit. 

This point of view was disputed by Glatzel (199^), who 
stressed that the fMimit is an artefact based on the disre- 
gard of gravity darkening and on the assumption of a uni- 
form brightness over the surface of rotating stars. Glatzel 
concludes that the Eddington factor has no effect on the 
critical rotation. This problem needs to be further exam- 
ined and this is what is done here. 

Another important issue is the dependence of the mass 
loss rates M on rotation velocities v (cf. Maeder and 
Me ynet |2000[ ). On one side, Lan ger (|1997| , |1998| ), Heger et 
al. fl2000| ), Meynet and Maeder ( pOPOf are using a relation 
M vs. v from Friend and Abbott (1986) which formally 



leads to infinite mass loss rates at break-up velocities. On 
the other side, Owocki et al. (1996], 1998), Owocki and 
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Gayley ( |l997| ) and Glatzel ( |1998|) show that even at ex- 



treme rotation the mass loss rates remain finite and do 
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not differ too much from the case of zero-rotation. There 
also, some further analysis is needed. 

In a rotating star, the total gravity is the sum of the 
gravitational, centrifugal and radiative accelerations: 



Stot — Seff + fifrad — Sgrav + fifrot + <?rad • 



(1.1) 



For purpose of clarity, we adopt the following definitions: 

- We speak of the Eddington or T-limit, when rotation ef- 
fects can be neglected and </rad + <7grav = 0, which implies 
that 



r = 



kL 



iircGM 



(1.2) 



In that case L = I-edd = AircGM/ k. The opacity k con- 
sidered here is the total opacity, unless we specify it dif- 
ferently (cf. Sect. 4.2). 

-The break-up or fi-limit is reached, for a star with an 
angular velocity Q at the surface, when the effective grav- 
ity g e f[ = 9grav + <7rot = and in addition when radiation 
pressure effects can be neglected. 

— The QT-limit is reached when the total gravity gtot = 0, 
with significant effects of both rotation and radiation. This 
is the general case, that we study here. It should lead to 
the two above cases in their respective limits. 

In Sect. 2, we examine the surface gravity, the Ed- 
dington factors and the limiting luminosity. In Sect. 3, 
the expression of the break-up velocities are considered, 
while in Sect. 4 we examine the mass loss rates. The equa- 
tion of the surface with account of rotation and radiative 
acceleration is discussed in an Appendix. 

2. Surface gravity, Eddington factors and limiting 
luminosity 

2.1. The von Zeipel Theorem close to the SIT -limit 



The von Zeipel theorem (1924) expresses that the radia- 
tive flux F at some colatitude i? in a rotating star is pro- 
portional to the local effective gravity g e ff. In a previous 



work (Maeder 1999), we have generalized this theorem to 
the case of shellular rotation proposed by Zahn ( 1992[ ). 
Shellular rotation results from strong horizontal turbu- 
lence which reduces the latitudinal dependence of rotation 
and makes the angular velocity constant on an isobar. 
Here, we shall consider the case of stars with shellular 
rotation, where the fir-limit may play a role. 

As shown by Langer ( |1997| ), stars close to the Edding- 
ton limit tend to develop convection in the outer layers (cf. 
also Maeder 1980). However, in the outer layers the con- 



vective flux is generally negligible and the main transport 
mechanism is rad iative transfer, a point also emphasized 
by Glatzel ( |l998|) . As a nu merica l example, in a 60 M 
model (Meynet and Maeder |200C|) at the end of the MS 



phase with log L/L Q = 5.89 and log T G ff =4.34 , the con- 
vective flux is negligible down to a fractional radius r/R 
of 0.85. In a 120 M Q with log L/L Q = 6.32 and log T cff 
=4.35 at the end of the MS, the convective flux is negligi- 
ble down to r/R = 0.63. Thus, the basic condition to apply 
the von Zeipel theorem to stars close to the fir-limit is 
fulfilled, since the flux is essentially radiative. 



The expression of the flux F (Maeder 199£) for a star 



with angular velocity f2 on the isobaric stellar surface (cf. 
Appendix) is 



F = - d^ flrf[1+CW] with 

M *=K i_ ^y and 



1 



f) e 



h t de 

S dr 



P 2 (costf). 



(2.3) 
(2.4) 
(2.5) 



There, p m is the internal average density, x — 
AacT 3 1 (inp) and \t is the partial derivative with respect 
to T . The quantity is defined by = 4, i.e. the ratio of 
the horizontal density fluctuation to the average density 
on the isobar, which is given by £ = -I^^P- where q is 



the average gravity on an isobar (cf. Zahn 1992). One has 
the thermodynamic coefficients 6 — — (d In p/d In T)p )M , 
Ht is the temperature scale height. The term which 
expresses the deviations of the von Zeipel theorem due to 
the bar oclini city of the star, is generally very small (cf. 
Maeder |l999| ). 

Let us emphasize that the flux is proportional to g e ff 
and not to gtot- This results from the fact that the equa- 
tion of hydrostatic equilibrium is ^f- = g e ff- The effect 
of radiation pressure is already counted in the expression 
of P, which is the total pressure. We may call M* the ef- 
fective mass, i.e. the mass reduced by the centrifugal force. 
This is the complete form of the von Zeipel theorem in a 
differentially rotating star with shellular rotation, whether 
or not one is close to the Eddington limit. 

2.2. Expressions of the gravity and of the local Eddington 
factor 

Let us express the total gravity at some colatitude d, tak- 
ing into account the radiative acceleration 



l„ n n($)F 

flrad = - V Pi ad = , 

P c 

thus one has with Eq. (1.1), (2.3) and (2.4) 



(2.6) 



9tot — g e f£ 



l - 



K (t?)£(p)[i + cm 

incGM (1 - 



n 2 



2ttGp„ 



(2.7) 



The rotation effects appear both in g e ^ and in the term in 
brackets. When we write k($), we mean that in a rotating 
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star, the local T e g and gravity vary with latitude and so 
does the opacity. We may also consider the local limiting 
flux. The condition g to t = in Eq. (1.1) with Eq. (2.6) 
for </rad allows us to define a limiting flux, 



fllm(0) 



re(tf) 



(2.8) 



From that we may define the ratio Ta(i3) of the actual 
flux F($) to the limiting local flux in a rotating star, 



r n (0) 



F(d) k(0)L(P)[1 + C(0)] 



AircGM 



(l- ° 2 

\/ 2TTG Pin ) 



(2.9) 



As a matter of fact, rn($) is the local Eddington ratio. 
For zero rotation r^f?) = T as given by Eq. (1.2). Using 
relation (2.9), we may write the Eq. (2.7) for the total 
gravity as 



gtot = <7eff [i-r n (0)] 



(2.10) 



This shows that the expression for the total acceleration 
in a rotating star is similar to the usual one, except that 
T is replaced by the local value Tn(^). Indeed, contrar- 
ily to expressions such as g to t = S^ff (1 — T) often found 
in literature, we see that the appropriate Eddington fac- 
tor (2.9) also depends on the angular velocity fi on the 
isobaric surface. 

From (2.9), we note that over the surface of a rotating 
star, which has a varying gravity and T e ff, Tn($) is the 
highest at the latitude where «;(#) is the largest, (if we 
neglect the effects of C($)j which is justified in general). 
If the opacity increases with decreasing T as in hot stars, 
the opacity is the highest at the equator and there the 
limit rn($) = 1 may be reached first. Thus, it is to be 
stressed that if the limit rn($) = 1 happens to be met at 
the equator, it is not because g e f? is the lowest there, but 
because the opacity is the highest ! Indeed, both depen- 
dences in g e ff have cancelled each other in the ratio given 
byEq. (2.9). 

2.3. The luminosity at the SIT -limit 

The fir-limit is reached, when the local Eddington ratio 
FnW = 1 at some colatitude The condition Tq^) = 1 
allows us to define a limiting luminosity Lor at the fir- 
limit, i.e. when both the effects of radiative acceleration 
and rotation are important. From (2.9) we have 



L 



AttcGM 



sir 



[i + co?)] 



i - 



fi 2 



2nGp n 



(2.11) 



It means that for a certain angular velocity fi on the iso- 
baric surface, the maximum permitted luminosity of a star 
is reduced by rotation, with respect to the usual Edding- 
ton limit (cf. S ect. 1). This conclusion was also reached by 
Glatzel (1998). In the above relation, k($) is the largest 



value of the opacity on the surface of the rotating star. For 
O-type stars with photospheric opacities dominated by 
electron scattering, the opacity k is the same everywhere 
on the star. For the equation of the surface discussed in the 
Appendix, the maximum value of = 0.361, (with 

more digits it is 0.360747). 

3. The break up velocities 

We have seen above that rotation may be considered as 
reducing the maximum possible luminosity for a star. An 
alternative way to consider the problem is to ask the 
question: what happens to the break-up velocity for a 
star close to the Eddington limit ? Most authors (Langer 
1994 ElU |1999| ; Lamers et al. |1999| Heger et al. |2000D 
write this critical velocity v CI n like 



GM 
~R~ 



(i-r) 



(3.12) 



This relation is true if we assume that the brightness of the 
rotating star is uniform over its surface, which is in contra- 
diction with von Zeipel's theorem. Surprisingly, some au- 
thors use this relation simultaneously with the von Zeipel 
theorem. Eq. (3.12), which we do not support, in agree- 
ment with Glatzel (1998), implies that the break-up ve- 
locity is reduced by the proximity to the Eddington limit. 
The problem needs to be further studied carefully. 

The critical velocity is reached when somewhere on the 
star one has gtot = 0, i.e. according to (2.10) 



9eff [i-r n (tf)] =0. 



(3.13) 



This equation has two roots. The first one f cr it,i is given 
by the usual condition g e ff = 0, which implies the equality 
Q 2 Rl h /(GM ) = 1 at the equator (cf. Eq. A2 in Appendix). 
This corresponds to an equatorial critical velocity 



^crit.l — fi Rcb 



2 GM 

3 R p b 



(3.14) 



R c h and R p h are respectively the equatorial and polar ra- 
dius at the break-up velocity and they obey to the surface 
equation. We notice that the critical velocity v cr it,i is in- 
dependent on the Eddington factor. To this extent, this 
is in agreement with Glatzel ( |1998| ). The basic physical 
reason for this independence is quite clear: the radiative 
flux decreases at the equator, when the effective gravity 
decreases. 

Equation (3.13) has a second root, which is given by 
the condition rn(i?) = 1. As seen above, this condition 
will in general be met at the equator first. We thus have 
to search for the corresponding critical velocity w cr it,2 for 
a given value of the stellar luminosity. This second root 
has to be compared to the first one. For given values of 
M and L, the lowest of the two roots u C rit,i an d v cr it,2 i s 
the significant one, since as soon as it will be reached the 
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matter at the surface of the star is no longer bound. The 
condition fnO?) = 1 gives, if we neglect CW> 



0.4 



k(<&)L(P) 
AircGM 

Let us write 



= 1 - 



n 2 



2wGp„ 



fi2 16 2t „, v 

2^ = 81* F (W) Wlth 



YM 



and 



2 _ 



(3.15) 



(3.16) 



(3.17) 



The quantity u> is the fraction of the angular velocity at 
the classical break-up given by Eq. (3.14). The density 
pmiu)) = M/V(uj), where the stellar volume V(uj) depends 
on rotation. The quantity V'iui) is the ratio of the actual 
volume of a star with rotation u> to the volume of a sphere 
of radius i? p b- V[w) is obtained by the integration of the 
solutions of the surface equation (A6) for a given value of 
the parameter u>. At break-up velocity u cr it,i , the value of 
V'{ijj) — 1.829, which gives the maximum value of 2 Jg p 
= 0.361. If we call r max the maximum Eddington ratio 
k({})L(P) I (AttcGM) over the surface (in general at the 
equator), Eq. (3.15) can thus be written 



- LO 2 V'(UJ) 

81 v ' 



= i-r r 



(3.18) 



For a given value of r max , one must search the value of 
u which satisfies this equation. This is easily obtained by 
solving numerically the surface equation. For a given large 
enough r max (i.e. larger than 0.639), the obtained w-value 
is lower than 1, and this implies a corresponding critical 
velocity Wcrit,2 given by 



81 1 



J crit 



a = nXM = 

9 



GM 



1 - 



i? 3 

^cb 



y crit,l 



(3.19) 



where R c (w) is the equatorial radius for a given value of 
the rotation parameter u>. Fig. [I] illustrates the results. We 
notice that the second root u C rit.2, expressed as a fraction 
of the first root v C rit,i given by (3.14), tends towards zero 
when the Eddington factor T max tends towards 1. The 
first root defined by (3.14) can still be written, but the 
second root defined by the condition r^(i?) = 1 is met 
first. This reduction of the critical velocity with respect to 

the classical expression only occurs for Eddington factors 

o 2 

larger than 0.639, since the maximum value of 2 -kGp ls 
0.361 (cf. Eq. 3.15). 

This second root results physically from both effects 
of rotation and radiation: for r max > 0.639, a zero value 
of 9tot ca n be achieved for non-extreme rotations. This 




0.4 O.i 

Fig. 1. The critical velocity u cr it,2 expressed as a fraction 
of u C rit,i plotted as a function of the Eddington factor 
r m ax, which is the largest value of the Eddington factor 
reached where the opacity is the largest over the stellar 
surface. We notice that when the Eddington factor tends 
towards unity, the critical velocity goes down to zero. 



enters through the reduction due to rotation of the effec- 
tive mass M*, which is the significant mass in the local 
Eddington factor To(^). 

For T max < 0.639, Eq. (3.18) has no solution. Physi- 
cally this means that when the star is sufficiently far from 
the Eddington limit, the reduction of the effective mass 
M* by rotation is not sufficient to bring rn(#) to 1. In 
that case, Eq. (3.13) has only one root given by the clas- 
sical Eq. (3.14). 

We h ope th at these resu lts cla rify the debate between 
Glatzel ( 1998 ) and Langer(1997). On one hand, we see 
that the claim by Langer that stars close to the Eddington 
limit have a lower rotation limit is correct, even if the Eq. 
(3.12) by Langer is not the right one. On the other side, 
Glatzel has claimed that the Eddington factor does not 
affect the break-up velocity, we see that this is true in 
general for most stars for the reasons given above, however 



for F max > 0.639 this statement does not apply. 

The moment when stars reach their critical velocities 
is far from being an academic one, since when this occurs 
large mass loss enhancements may result, a point which is 
examined below. 



4. The mass loss rates as a function of ft and T 
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4--1- Present context 

The effects of rotation on the mass loss rates have been 
studied both observationally and theoretically. Observa- 
tionally, very large changes of the M-rates, i.e. up to 2- 3 
orders of a magnitude, were suggested by Vardy a ( 1985p . 
However, Nieuwenhuijzen and de Jager ( 1988j ) claimed 
with reason that the correlation found by Vardya was 
largely the reflect of the distribution of the mass loss 
rates and rotational velocities over the HR diagram. When 
disentangling the various effects, Nieuwenhuijzen and de 
Jager found much smaller effects of rotation. However, 
they noticed that the M-rates of the Be-stars are larger 
by about a factor 10 2 . Since Be-stars are fast rotating 
stars, we may wonder whether the effects of rotation on 
the mass loss rates are really so negligible, as these last 
authors considered them. Certainly further observational 
studies are also required. 



On the theoretical side, Pauldrach et al. (1986), Friend 



and Abbott ( 1986 ) find only a moderate increase of the 
M-rates, of about 30 % for v = 350 km/s. Friend and 
Abbott find an increase of the M-rates which can be fitted 



by the relation (Langer 1998 ; Heger and Langer 1998 ) 
M{v) = M(v = 0) (7^^) ( 4 - 2 °) 



with £ = 0.43. This expression, often used in evolutionary 
models, is based on wind models which do not account for 
the von Zeipel theorem. We notice that Eq. (4 .20) diverges 
at break-up, while as shown by Glatzel (1998) and Owocki 
et al. (1996, 1998), the stellar mass loss rates should not 
diverge at the O-limit, (however see Sect. 4.3). 

When a proper account of the gravity darkening is 
made, there are two main terms contributing to the 
anisotropic mass loss rates from a rotating star (cf. Maeder 
1999| ). 1) The "g e gf-effect" which favours polar ejection, 



since the polar caps of a rotating star are hotter. 2) The 
"opacity or K-effect" , which may favour an equatorial ejec- 
tion, when the opacity is large enough at the equator due 
to the lower T e g. 

In O-type stars, since opacity is due mainly to the 
T-independent electron scattering, the g e ff-effect is likely 
to dominate, raising a fast highly ionized polar wind. In 
B- and later type stars, the opacity effect may favour a 
dense equatorial wind and ring formation, with low ter- 
minal velocities and low ionization. Recently Petrenz and 
Puis ( 200Cj ) have constructed 2-D models of line driven 
winds for rotating O-type stars. They found that the mass 
loss from hot star is essentially polar due to the g e ff-effect. 
Quantitatively their results differ very little from previous 
works by Pauldrach et al. ( |1986[ ). 

4-2. Expression of the M-rates as a function ofQ andT 

It is worth to further examine the consequences of the 
above results on the dependence of the mass loss rates 



on rotation. According to the radiative wind theory (cf. 
Castor et al. 1975; Pauldrach et al. 1986: Kudritzki et al. 



|1989| Puis et al. 1996), we may write the mass loss fluxes 



AM/ Act by surface elements Act 



AM(tf) 
Act 



(ka 



1/q 



F^) 1/a 9 ± tat "W (4-21) 



where k and a are the force multiplier parameters. At 
some temperatures, the ionisation equilibrium of the stel- 
lar wind is changing abruptly and so does the opacity of 
the plasma. Consequently, the values of the force multi- 
pliers undergo rapid transitions for certain values of T e g, 
particularly at 21 000 K and maybe also at 10 000 K (cf. 
Lamers et al. 1995; Lamers 1997). Such fast transitions 



of the wind properties are called by Lamers a bi-stability 
of the stellar winds, since near the transition limit the 
wind can exist in two states. There are both empirical 
and theoretical determinations of a, however they lead to 



rather different values (cf. Lamers et al. 1995). The em- 
pirical ones, based on the values of the observed terminal 
velocities, are in general smaller than the theoretical es- 



timates. As empirical values, Lamers et al. (1995) obtain, 
for example, a = 0.52 for 4.70 > logT eff > 4.35, (type 
B1.5 or earlier); a = 0.24, 0.21, 0.17, 0.15 for logT off = 
4.30 (type B2.5), 4.20 (B5), 4.00 (B9.5), 3.90 (A7) respec- 
tively. These transitions may produce jumps in the mass 
loss rates, with the high rates on the low side of the tran- 
sition. As T e ff is decreasing from the pole to the equator, 
one may thus expect, on the surface of a fast rotating star 
of type B or later, the occurence of some bi-stability limits 
and the corresponding variations of the force multipliers 
and of the mass loss rates. In a star, where a bi-stability 
limit is crossed at some latitude, a steep increase of the 
mass flux will happen between this latitude and the equa- 
tor, possibly leading to a huge equatorial ring. 

With the expressions of the flux (2.3) and of gtot 
(2.10), we get for the mass flux 



AM(tf) 
Act 



.4 



L(P) 



4ttGM*(P) 
with A = (ka) 



geff[l + Cg]f 

(l-roW)^ 1 

1 — Of" 



(4.22) 



with M* given by Eq. (2.4). We notice the gravity-effect, 
which favours mass loss at the pole, where the total gravity 
is higher, and the K-effect which favours high mass loss 
where a is small. The proximity to the Eddington limit 
will enhance the mass flux due to the term Fq($), while 
rotation enhances the mass flux through both the terms 
M* and Tn(^). In the theory of radiatively driven winds, 
the total opacity at a given optical depth is expressed with 
the force multipliers in terms of the electron scattering 
opacity k cs . This means that in Eq. (4.22), rn($) is just 
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r n (0) 



l{p)[ i + cm 

AitcGM ( 1 



2ttGp„ 



(4.23) 



The dependence on latitude of Tq(-&) would only come 
through the term 

4-3. Dependence of the global mass loss rates on rotation 

Let us estimate how the global mass loss rates depend 
on the rotation velocities and on the proximity of the £IT 
limit. For that we henceforth neglect the small corrective 
term ((d) in the expression of the flux. We have, if £(w) 
is the total surface 



M 
EH 



.4 



L(P) 



9eS 



(i-r n )t 



.4 



L(P) 



ffcir 



(1 



(4.24) 



since Tq is independent on and with appropriate a- 
and ^-values. For the average effective gravity, we have 



ffcff 



/ / geff • da- _ AitGM* 



EH 



(4.25) 



after integration over the stellar surface which is an isobar 
(cf. Appendix). This leads to the following expression for 
the total mass loss rate from the star 



M 



AL(py. 



(AttGM)' 



n 2 



2ttGp„ 



(i-r„)*- 



(4.26) 



This relation expresses how the total mass loss rate from 
a star depends on mass, luminosity, Eddington factor and 
rotation, (see Fig. Al for a simple expression of the rota- 
tion parameter). If we omit rotation, Eq. (4.26) is identical 
to the typical relations used in literature (cf. Pauldrach ct 
al. 1986| ; Lamcrs [l997 ). The amplitude of the effects very 
much depends on the opacity and in particular on the 
value of the force multiplier a. 

Let us consider a rotating star with angular velocity f2 
and a non-rotating star of the same mass M at about the 
same location in the HR diagram. The ratio of their mass 
loss rates can be written, 



M(0) 
M(0) 



(i-r)< 



i 



n 2 



2irGp„ 



(4.27) 



(l-r„); 



where T is the Eddington ratio corresponding to electron 
scattering opacity for the non-rotating star. From Eq. 
(4.23), we have the relation 



ro = 



n 2 

2wGp n 



(4.28) 



Table 1. Values of T at the end of the MS for various 
initial stellar masses, and of the ratios M(£Y)/M(0) of the 
mass loss rates for a star at break-up rotation to that 
of a non-rotating star of the same mass and luminosity 
at logTcff > 4.35, at logT off = 4.30, 4.00 and 3.90 The 
empirical force multipliers a by Lamers et al. (1995) are 
used. 



M ini 


r 


Mill) 

M(0) 

a = 0.52 


il'/lS!) 

M(0) 

a = 0.24 


M(0) 

a = 0.17 


M(0) 

a = 0.15 


120 


0.903 


oo 


oo 


oo 


oo 


85 


0.691 


oo 


oo 


oo 


oo 


60 


0.527 


3.78 


96.2 


1130 


3526 


40 


0.356 


2.14 


13.6 


55.3 


106.0 


25 


0.214 


1.76 


7.02 


20.1 


32.6 


20 


0.156 


1.67 


5.87 


15.2 


23.6 


15 


0.097 


1.60 


5.04 


12.1 


18.1 


12 


0.063 


1.57 


4.68 


10.8 


15.8 


9 


0.034 


1.54 


4.41 


9.8 


14.2 



We get finally 

M(O) _ (1-r)"" 1 
M(0) ~" 



1 - 



n 2 

2irGp„ 



-r 



(4.29) 



If £1 = 0, this ratio is of course equal to 1. This ratio, which 
is the main result of this work, can also be expressed with 
the ratio w/fcrit,i of the rotational velocity v to the critical 
velocity given by the usual Eq. (3.14). From Eq. (A7) in 



the Appendix, we have 



2ttGp„ 



9 v 2 



over a large range 



of values (cf. Fig. Al), thus one can write 



M(0) 



(1 



1 



4 (_v_\2 

9 Writ, 1 I 



(4.30) 



For a star with a small Eddington factor, it simplifies to 
M(Q) 1 



M(0) 



1 



V l>crit,l ' 



(4.31) 



Equation (4.31) shows that the effects of rotation on 
the M-rates remain moderate in general. This is in agree- 
ment with the results by Owocki et al. (|1996 ), by Owocki 
and Gayley ( |1997| ) and by Glatzel ( |1998| ), and also with 
more elaborate non-LTE 1-D and 2-D models by Paul- 
drach et al. ( |1986| ), Petrenz and Puis ( |2000| ). However, this 
is only true for stars far enough from the Eddington limit. 
When r is significant, rotation may drastically increase 
the mass loss rates as shown by (4-29) or (4-30). This is 
particularly the case for low values of a, i.e. for stars with 
log T c ff < 4.30. In the extreme cases where T > 0.639, 
a moderate rotation may even make the denominator of 
(4.29) or (4.30) to vanish, thus leading to extreme mass 
loss. 



A. Maeder and G. Meynet: Stellar Evolution with Rotation 



7 



Table 1 shows some numerical results based on Eq. 

(4.29) . For different initial stellar masses in the Geneva 
models at Z =0.02 (Schaller et al. |l992j ), the T factors at 
the end of the Main Sequence (MS) phase are given, as well 
as the predicted ratios M(f2)/M(0) of the mass loss rates 
for a star at break-up rotation and for a non-rotating star 
of the same mass and luminosity. These ratios are given at 
logT cff > 4.35 (a = 0.52), at logT eff = 4.30 (a = 0.24), at 
logT off = 4.00 (a = 0.17) and at logT cff = 3.90 (a = 0.15) 
for the same value of V. This covers the range of the typical 
T e g of OB and LBV stars, the differences with T c ff result 
from the differences in the a-parameter. The indication 
oo in Table 1 means that the bracket term in (4.29) or 

(4.30) may vanish at maximum rotation, which leads to 
extreme mass outflows. 

The ratios M(f2)/Af(0) keep quite moderate even at 
extreme rotation for MS stars up to 40 Mq, while for MS 
stars above 60 M they can become very large. These ra- 
tios may also be very large for B-type supergiants and 
LBV stars. In particular, we notice that for stars close 
to the Humphreys-Davidson limit the ratios M(f2)/A/(0) 
may diverge. Such stars are typically at the fir-limit. For 
log T e g < 4.30, the force multiplier a is also very small, 
which favours extreme mass loss. On the whole, it is strik- 
ing that the domain where M (fi) j M (0) has the possibility 
to diverge so closely corresponds to the observed domain 
of LBV stars. The present results will enable us to better 
specify the changes of the M rates in massive star models. 



imity to the Eddington limit. As shown by Table 1, both 
effects produce steep enhancements of the mass loss rates, 
especially for lower T B g since a is lower. This may explain 
the very large mass loss rates for LBV stars, blue and 
yellow supergiants (cf. de Jager et al. 1988). If the ratio 



r = 



4vcGM 

reached for reduced rotation velocities 



is bigger than 0.639, the break-up limit is 

as illustrated by 
Fig. 1. Then, extremely high mass loss rates may occur, a 
situation likely corresponding to the case of the LBV stars 
and maybe also to some WR stars. 

Some words of caution are necessary. The M-rates 
given here are the values predicted in the framework of the 
radiative wind theory. It is probable that close to break-up 
several other effects not included here may intervene, such 
as important horizontal fluxes, formally vanishing T— and 
P- gradients, instabilities, etc... Also, we may point out 
that if the flux vanishes, the radiative wind theory should 
not apply. Thus, for the detailed physics of the break-up, 
more complex analyses are certainly needed. 

Finally, we note that it was generally believed that in 
addition to L and T e g, the mass loss rates only depend 
on metallicity Z. We see here another dependence which 
is quite significant and may introduce some scatter in the 
values of the M-rates. Thus, we may expect that for a 
given initial mass the evolution is very different accord- 
ing to rotation, due to both rotational mixing, meridional 
circulation (Maeder and Zahn 1998 ) and to the induced 
differences in the mass loss rates. 



5. Conclusion 

We conclude that the concept of an fir-limit reached dur- 
ing the evolution of the most massive stars is not an arte- 
fact, but the existence of this limit is confirmed by consis- 
tent developments based on the von Zeipel theorem. How- 
ever, we emphasize that the expression currently used for 
the critical velocity is not correct. We have also clarified 
the dependence of the mass loss rates on the rotation ve- 
locities in the general case. 

We can make the following remarks on the various lim- 
its: 

—1. The r limit : The mass loss rates grow steeply 
as the Eddington limit is approached, even in absence of 
rotation. This is a well known result of the classical wind 
theory. 

2. The fi— limit : We see that the case of only rota- 
tional effects does not apply for O-type stars and even for 
the early B-type stars, since they always have a signifi- 
cant T-value. Only for spectral types later than B3 on the 
MS, the F term can be ignored. In the framework of the 
radiative wind theory, the growth of the mass loss rates 
remains limited. 

—3. The fir— limit : This general case is met for ro- 
tating OB stars, LBV stars, supergiants and Wolf-Rayet 
stars, because both T and rotation are important. The 
bracket in (4.29) is reduced by rotation and by the prox- 



Appendix A: The equation of the stellar surface 
in a rotating star with high 
radiation pressure 

Shellular rotation, with an angular velocit y fi co nstant on 
horizontal surfaces, was proposed by Zahn ( 1992 ). This ro- 
tation law results from strong horizontal geostrophic-like 
turbulence which homogeneizes rotation on th e hor izon- 
tal surfaces. As noted by Meynet and Maeder ( 1997 ), the 
isobars for shellular rotation are identical to the equipo- 
tentials of the conservative case, which are 



GM 1 

* = — -rr + -fiV(tf) sin 2 = const 
r(v) 2 

The components of the effective gravity are 



<9$ 2 
g e ft,r = — + fi r sin $ 
or 



19$ , 
ffcff i? = - ttt r + fi t sin $ cos i9 
r ov 



(A.l) 



(A.2) 



where $ = GM/r. In vectorial form, one can write 

VP = -pg eff = -p(V* - r 2 sin 2 & fiVfi) . (A.3) 
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Fig. A.l. The rotation parameter 



as a function 



2itGp„_ 

of the square of the ratio of the rotational velocity to the 
critical velocity (3.14). The maximum value of 2 ^ l Gp = 
0.361, as seen above. 



This equation is interesting: it shows that if is constant 
on isobars, * is also constant on isobars. Moreover, for a 
motion ds, the equation of the surface must satisfy 



9tot • ds = . 
Since we have 
<7tot = Seff [1 - ToJ 



(A.4) 



(A.5) 



this implies that the surface is perpendicular to gtot , thus 
to g e ff and to the P — gradient. The surface is an isobar, 
also if the radiation pressure is important. The equation 
of the surface can therefore be represented by Eq. (Al) 
and the procedure to calculate is quite conventional, i.e. 



1 i ^ 2 2 • 2 a 
1 U X sin V 

x 27 



(A.6) 



with x — At break-up, the equatorial radius R e \y 
equals 1.5 times the polar radius R p b- 

In the above demonstration, there is no need of the as- 
sumption £1 = f2(r), we just need the assumption that f2 
is constant on horizontal surface, which is less restrictive. 
This is true, whether the radiative acceleration is impor- 
tant or not. 

The critical rotation parameter naturally appearing in 



this work was the ratio 



Account must be given to 
the change of the average density of the star with rotation. 



2ttGp„ 



We can express this ratio by (3.16) or in term of the ac- 
tual rotational velocity v and of the critical velocity w cr it i 
(3.14), 



ti 2 



9 ^crit.l 



(A.7) 



The relation between 



2ttGp„ 



and the ratio 



T) 2 



— is illus- 

"crit.l 

has a limited 



trated in Fig. Al. The product V'(u) 
range of variation, being equal to 1 for zero rotation and 
to 0.813 at break-up velocity. This means that for a crude 
estimate at low or moderate velocities, one may just ignore 
this product in Eq. (A7). 
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